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Abstract. KASCADE-Grande is a multi-detector 
experiment at Forschungszentrum Karlsruhe, Ger- 
many for measuring extensive air showers in the 
primary energy range of 100 TeV to 1 EeV. This 
contribution attempts to provide a synopsis of the 
current results of the experiment. In particular, the 
reconstruction of the all-particle energy spectrum in 
the range of 10 PeV to 1 EeV based on four different 
methods with partly different sources of systematic 
uncertainties is presented. Since the calibration of the 
observables in terms of the primary energy and mass 
depends on Monte Carlo simulations, we compare the 
results of various methods applied to the same sample 
of measured data. In addition, first investigations on 
the elemental composition of the cosmic particles 
as well as on tests of hadronic interaction models 
underlying the analyses are discussed. 

Keywords: High-energy cosmic rays, extensive air 
showers, KASCADE-Grande 

I. Introduction 

The main goal of experimental cosmic ray research is 
the measurement of the primary energy spectrum and the 
elemental composition. That comprises important hints 
to understand the origin, acceleration and propagation 
of energetic cosmic particles. This task can be done 



directly or indirectly, depending on the energy of the 
primary particle. At high energies, above 10 15 eV the 
energy spectrum must be determined indirectly from the 
measured properties of extensive air showers (EAS) that 
cosmic rays induce in the Earth's atmosphere. 

Depending on the experimental apparatus and the 
detection technique, different sets of EAS observables 
are available to estimate the energy of the primary 
cosmic ray Q. In ground arrays the total number of 
charged particles in the shower and the corresponding 
density at observation level are commonly employed Q. 
However, the muon content of EAS can play also an 
important role. One reason in favor of this observable 
is that in an air shower muons undergo fewer atmo- 
spheric interactions than the electromagnetic or hadronic 
components (dominant component for vertical EAS) and 
exhibit in consequence less fluctuation compared to the 
more abundant electromagnetic component. 

In KASCADE-Grande both components, the muon 
and the electromagnetic are measured and both, together 
with their correlation on a single-event-basis are used to 
derive the energy spectrum and composition of cosmic 
rays in the range from 10 16 to 10 18 eV. 
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Fig. 1: Layout of the KAS CADE- Grande experiment: 
The original KASCADE, the distribution of the 37 
stations of the Grande array, and the small Piccolo 
cluster for fast trigger purposes are shown. The outer 
12 clusters of the KASCADE array consist of /i- and 
e / 7-detectors, the inner 4 clusters of ej 7-detectors only. 



II. KASCADE-Grande 

Main parts of the experiment are the Grande array 
spread over an area of 700 x 700 m 2 , the original KAS- 
CADE array covering 200 x 200 m 2 with unshielded and 
shielded detectors, and additional muon tracking devices. 
This multi-detector system allows us to investigate the 
energy spectrum, composition, and anisotropies of cos- 
mic rays in the energy range up to 1 EeV. The estimation 
of energy and mass of the primary particles is based on 
the combined investigation of the charged particle, the 
electron, and the muon components measured by the 
detector arrays of Grande and KASCADE. 

The multi-detector experiment KASCADE EJ (lo- 
cated at 49.TN, 8.4°E, HOma.s.l.) was extended to 
KASCADE-Grande in 2003 by installing a large array 
of 37 stations consisting of 10 m 2 scintillation detectors 
each (fig. Q]). KASCADE-Grande H provides a sensi- 
tive area of « 0.4 km 2 and operates jointly with the 
existing KASCADE detectors. The joint measurements 
with the KASCADE muon tracking devices are ensured 
by an additional cluster (Piccolo) close to the center of 
KASCADE-Grande for fast trigger purposes. 

While the Grande detectors are sensitive to charged 
particles, the KASCADE array detectors measure the 
electromagnetic component and the muonic component 
separately. The muon detectors enable to reconstruct 
the total number of muons on an event-by-event basis 
also for Grande triggered events. The Muon Tracking 
Detector (MTD) 01 registers muons above an energy 
threshold of 800 Me V. At the MTD the directions of 
muon tracks in EAS are measured with an excellent 
angular resolution of « 0.35°. These directional data 
allow to investigate the longitudinal development of 
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Fig. 2: Example of a typical KASCADE-Grande event. 
Displayed are the measured densities per station for 
charged particles of the Grande array and averaged over 
stations in 20 m rings for the KASCADE muon array, 
respectively and the corresponding fits of the lateral 
distribution functions. 



the muonic component in showers which is used as a 
signature of the development of the hadronic EAS core. 

III. Reconstruction 

Basic shower observables like the core position, angle- 
of-incidence, and total number of charged particles are 
provided by the measurements of the Grande stations. 
By means of Monte Carlo simulations a core position 
resolution of ~ 5 m, a direction resolution of « 0.7°, and 
a resolution of the total particle number in the showers of 
w 15% is expected |6]. The total number of muons (N^ 
resolution w 25%) is calculated using the core position 
determined by the Grande array and the muon densities 
measured by the KASCADE muon array detectors Q. 
Fig. [2] shows a typical event of KASCADE-Grande 
and its reconstruction; i.e. the functions fitted to the 
measured lateral particle density distributions. Full ef- 
ficiency for triggering and reconstruction of air- showers 
is reached at primary energy of « 2 • 10 16 eV, slightly 
depending on the cuts needed for the reconstruction of 
the different observables. 

For a subsample of events collected by the Grande 
array it is possible to compare on an event-by-event basis 
the two independent reconstructions of KASCADE and 
Grande. This provides the unique opportunity of evalu- 
ating the reconstruction accuracies of the Grande array 
by a direct comparison with an independent experiment. 
Since the KASCADE array is much more dense, the 
values of reconstructed observables from KASCADE 
are taken as reference. The subsample is obtained ac- 
cordingly to the following selection criteria: maximum 
energy deposit in one of the Grande stations close to 
KASCADE; core position within a circle of 90 m radius 
from the KASCADE center; zenith angle less than 40°. 
By means of such a comparison (see fig. the Grande 
reconstruction accuracies are found to be for the shower 
size: systematic uncertainty < 5%, statistical inaccuracy 
< 15%; for arrival direction: a w 0.6°; for the core 
position: a « 5 m. All of them are in good accordance 
with the resolutions obtained from simulations. 
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Fig. 3: Scatter plot of the shower sizes (charged parti- j s s hown. 
cles) reconstructed by KASCADE (x-axis) and Grande 
(y-axis). 



Fig. 4: Mean distribution of true energy vs. muon num- 
ber for iron nuclei and = 20.5° — 26.5° calculated with 
MC simulations. The fit with formula E[GeV] = /? • 



IV. All-particle energy spectrum 

A. Strategy 

Applying different methods to the same data sample 
has advantages in various aspects: one would expect 
the same result for the energy spectrum by all methods 
provided the measurements are accurate enough, when 
the reconstructions work without failures, and when the 
Monte Carlo simulations describe correctly the shower 
development. But, the fact that the various observables 
exhibit substantial differences in their composition sen- 
sitivity hampers a straightforward analysis. However, 
investigating results of different methods can be used 
to 

• cross-check the measurements by different sub- 
detectors; 

• cross-check the reconstruction procedures; 

• cross-check the influence of systematic uncertain- 
ties; 

• test the sensitivity of the observables to the elemen- 
tal composition; 

• test the validity of hadronic interaction models 
underlying the simulations. 

B. Methods 

The estimation of the all-particle energy spectrum is 
presently based on four different methods using different 
observables of KASCADE-Grande: 

• N c h -method: The reconstructed charged particle 
shower size per individual event is corrected for 
attenuation in the atmosphere by the constant in- 
tensity cut method and calibrated by Monte Carlo 
simulations under the assumption of a dependence 



Eo oc N^ h and a particular primary composi- 
tion n. 

• A^ -method: The reconstructed muon shower size 
per individual event is corrected for attenuation and 
calibrated by Monte Carlo simulations under the as- 
sumption of a dependence Eq oc (see example 
in fig. [4]) and a particular primary composition 0. 

• N c h — A/^-method: This method combines the infor- 
mation provided by the two observables. By help of 
Monte Carlo simulations a formula is obtained to 
calculate the primary energy per individual shower 
on basis of N c h and iV M . The formula takes into 
account the mass sensitivity in order to minimize 
the composition dependence. The attenuation is 
corrected for by deriving the formula for different 
zenith angle intervals independently and combining 
the energy spectrum afterwards iflQl . 

• ^(SOOj-method: The reconstructed particle density 
at the specific distance to the shower axis of 500 m 
per individual event is corrected for attenuation and 
calibrated by Monte Carlo simulations under the 
assumption of a dependence i^o oc 5'(500) Q!s ( 500 ) . 
The distance of 500 m is chosen to have a minimum 
influence of the primary particle type, therefore 
a smaller dependence on primary composition is 
expected ifTTl . 

In figure \5\ the resulting spectra are compiled. Owing 
to the different procedures, the results for the first two 
methods are shown under proton and iron assumption, 
respectively, only, whereas for the other two methods 
the resulting all-particle energy spectra are displayed. 
Figure [6| shows the same results but with the flux 
multiplied by a factor of E 3 - . 
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Fig. 5: Reconstructed all-particle energy spectrum by 
four different methods applied to KAS CADE- Grande 
data. Given are also the energy resolutions for the 
methods. 
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Fig. 6: Same as figure [2 but the flux multiplied by Eq. 
Values for the uncertainty in the flux determination are 
given for the different methods. 



C. Systematic uncertainties and attenuation 

The application of the different methods allows us 
to compare and cross-check the influence of various 
sources of systematic uncertainties. The N c h -method 
uses the basic measurements of the Grande array only, 
resulting in a high accuracy of N c h with better than 15% 
over the whole range of shower size, without any energy 
dependent bias. But, using only one observable, there is a 
large dependence on the primary elemental composition, 
reflected by the distance between the spectra obtained 
for proton and iron assumption at the calibration. The 
N^ -method on the other hand is based on the muon 
shower size, which can be derived less accurately (25% 
with a small bias dependent on the distance of the 
shower core to the muon detectors which is corrected 
for), but with much less composition dependence. The 
N c h-N^ -method, owing to the combination of the recon- 
struction uncertainty of two variables, shows basically 
a larger uncertainty in the reconstruction, but this is 
compensated by taking into account the correlation of 
these observables at individual events. Furthermore, by 
this procedure the composition dependence is strikingly 
decreased. The 5(500) value by construction yields a 
larger uncertainty of the variable reconstruction, but has 
also a minor composition dependence. 

For all methods, the energy resolution is estimated 
using full Monte Carlo simulations and comparing the 
reconstructed with the simulated primary energy (for 
instance figure [5] gives the numbers for an energy of 
Eq = 10 17 eV). Values of systematic uncertainties in the 
flux determination for the different methods are shown in 
fig. [6] (again for Eq = 10 17 eV). To a large amount these 
uncertainties are due to the reconstruction of the ob- 
servables, but there are additional sources of systematics 
which belong to all methods: e.g., concerning the Monte 
Carlo statistics, the assumed Monte Carlo spectral slope, 
or the fits of the calibration procedures. The different 



attenuation (and its handling to correct for) of the various 
observables (A(N ch ) « 495 ± 20g/cm 2 ; A(N^) w 
1100 ± 100g/cm 2 ; A(S(500)) « 347 ± 22g/cm 2 at 
Eq = 10 17 eV) however, lead again to slightly different 
contributions to the total systematic uncertainty. The 
total uncertainties (energy resolution and systematics) 
for the various methods are discussed in refs. O, (9), 
ifTOl . lITTIl and can be displayed as a band surrounding 
the reconstructed energy spectrum (e.g., see fig. 0. 

D. Discussion 

Taking into account the systematic uncertainties, there 
is a fair agreement between the all-particle energy spec- 
tra of the different applications (fig. 0). 

Of particular interest is the fact that by using N c h, 
the iron assumption predicts a higher flux than the 
proton assumption, whereas using the opposite is 
the case. That means that the 'true spectrum has to 
be a solution inside the range spanned by the two 
methods. If there is only the possibility of applying 
one method, then there is a large variance in possible 
solutions (everything within the range spanned by proton 
and iron line, not even parallel to these lines). However, 
more detailed investigations have shown, that a structure 
in the spectrum or a sudden change in composition 
would be retained in the resulting spectrum, even if the 
calibration is performed with an individual primary, only. 
Interestingly, over the whole energy range there is only 
little room for a solution satisfying both ranges, spanned 
by N c h and N^, and this solution has to be of relatively 
heavy composition - in the framework of the QGSJet- 
II hadronic interaction model. The narrower range for a 
solution provided by the A/^ -method compared to N c h 
confirms the finding of KASCADE that at sea-level the 
number of mostly low-energy muons N^ is a very good 
and composition insensitive energy estimator. 

The results of the composition independent N c h- 
N^-, and S (5 00) -methods lie inside the area spanned 
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Fig. 7: Lateral distribution of muons compared to the Fig 8: Muon density as a funct i n of the total number 
predictions of QGSJet-II and EPOS 1.61. of electrons compared to the predictions of QGSJet-II 

and EPOS 1.61. 



by the composition dependent methods, which is a 
very promising result. The £(5 00) -method results in 
a slightly higher flux than the N^-N^ -method, but 
the two spectra are consistent taking into account the 
systematic uncertainties. 

All the discussed results show a smooth all-particle 
energy spectrum without a clear hint to a distinct struc- 
ture over the whole energy range from 10 PeV to 1 
EeV. Another conclusion is that, taking into account the 
systematic uncertainties for all methods, the underlying 
hadronic interaction model (QGSJet-II/FLUKA) is in- 
trinsically consistent, i.e. the correlation between the dif- 
ferent observables, respectively the particle components 
can describe the global features of our measurements. 

V. Hadronic interaction models 

By now, for all the considerations the models QGSJet- 
II and FLUKA d, C3, d have been used, only. 
Other interaction models would probably change the 
interpretation of the data. In the following it is briefly 
described how we investigate the influence of hadronic 
interaction models on the interpretation of the measured 
data. 

A. Muon densities 

Local muon densities measured with the KASCADE 
array are studied in their correlation with distance from 
the shower axis and the total number of electrons in the 
shower lfT5l . 

Figure [7] shows the mean muon density as a function 
of the distance from the shower axis compared to the 
predictions of QGSJet-II and EPOS 1.61. Both hadronic 
interaction models include the data within the proton 
and iron limits for the entire range of distances from 
100 to 750 meters. Not only the absolute density but 
also the slope of the LDF is of interest for checking the 
validity of hadronic interaction models. Considering an 
equal probability trigger for protons and iron primaries 
as a function of distance from the shower axis, one 
should expect the measured LDF to be parallel to 



primaries of pure composition. Note that the LDFs of 
simulated proton and iron showers are parallel. However 
the measured LDF is neither parallel to the QGSJet-II 
nor to the EPOS 1.61 curves. It shows that the slope of 
the LDF cannot be well described by either model. 

Figure [8] shows the evolution of the mean muon 
density as a function of N e . The calculations done with 
QGSJet-II and EPOS 1.61 using proton and iron nuclei 
as primary particles embrace the data in the entire range 
of 5 < lg(N e ) < 8. 

Nevertheless, both figures [7] and [8] show that EPOS 
1.61 would require a very light primary composition in 
order to fit the data. On the other hand, QGSJet-II could 
fit the data with an intermediate primary abundance 
between proton and iron nuclei. 

B. All-particle energy spectrum 

To study effects on the all-particle energy spectrum 
we investigated the influence of the hadronic interac- 
tion model exemplary by performing the A^-method 
based on simulations with the hadronic interaction model 
EPOS vers. 1.61 d in addition to QGSJet-II. As the 
Monte Carlo statistics is limited in case of EPOS, 
both spectra were obtained by generating the calibra- 
tion curve with an equally mixed composition of five 
primaries (H,He,C,Si,Fe). Figure [9] compares the all- 
particle energy spectrum obtained with the KASCADE- 
Grande data set for both cases. The interpretation of the 
KASCADE-Grande data with EPOS leads to a signif- 
icantly higher flux compared to the QGSJet-II result. 
Though we know, that version 1.61 of the EPOS model 
is not fully consistent with air shower data (in particular, 
it cannot describe the correlation of hadronic observables 
with the muon and electron content of the EAS 1 17 ]) this 
example shows that applying and comparing various re- 
construction methods on the same data set will be useful 
for a better understanding of the interaction processes in 
the air shower development. 
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Fig. 9: Reconstructed all-particle energy spectrum with 
the N c h -method and the calibration function obtained by 
assuming mixed composition, but based on two different 
hadronic interaction models. 



C. Muon Tracking Detector 

More general investigations of the characteristics of 
the hadronic interaction models are provided by studies 
of the muon tracking data lfT8l . Main observables of 
interest are the mean muon radial angles and the mean 
muon pseudo-rapidities |[T9lL which are reconstructed 
in a given distance range from the shower core, and 
which are quantities sensitive to the development of 
the extensive air shower. Vertical showers (0 < 18°) 
with size lg7V e > 6 have been selected. The lateral 
distributions of those two quantities have been obtained 
and compared with the results reconstructed out of the 
simulated data for two primary species: proton and iron. 
The simulations were done with CORSIKA fl2ll code 
using the QGSJet-II model [ 13 1 for high energy inter- 
actions above 200 GeV and FLUKA2006 d below 
that energy. The comparisons are done in limited ranges 
of muon distances to the core, where the saturation 
effects (seen below 150 m) and trigger inefficiencies 
(seen above 400 m) are not present. For the mean radial 
angle (fig. [TO]) we can conclude that the experimental 
data are compatible with the CORSIKA simulations 
using the QGSJet-II - FLUKA model combination, as 
the data points are in between the simulated ones. The 
lines are linear fits to the simulation results. The error 
bars in the simulations are still very large and the 
number of simulated data will be increased, thus the 
results are marked as preliminary. We also notice that 
the experimental data tend to be positioned closer to the 
line of iron initiated showers rather than proton ones, like 
for the muon densities shown in fig. [7] In fig. [TTJ where 
lateral distributions of the mean muon pseudorapidity 
for the same data set are compared, one can also 
conclude that the experimental data points are bracketed 
by the simulated distributions, however, the data points 
here are closer to the proton simulation results. This 
inconsistency may be an indication of the features of the 
models. The mean radial angle distribution suggests that 
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Fig. 10: Reconstructed lateral distribution of the mean 
radial angle compared with CORSIKA simulation results 
for proton and iron primaries. Lines are fits to the 
simulations. 
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Fig. 11: Reconstructed lateral distribution of the mean 
muon pseudorapidity compared with CORSIKA simula- 
tion results for proton and iron primaries. Lines are fits 
to the simulations. 



the transverse momentum of pions produced in hadronic 
interactions is reproduced by the models, but fig. Qj] is 
a hint that the simulations provide relatively more high 
rapidity pions. 

VI. Composition 

The basic goal of the KAS CADE- Grande experiment 
is the determination of the chemical composition in 
the primary energy range 10 16 — 10 18 eV. Like for the 
reconstruction of the energy, again several methods using 
different observables will be applied to the registered 
data, in order to study systematic uncertainties. However, 
the influence of predictions of the hadronic interaction 
models has a much larger influence on the composi- 
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tion than compared to the primary energy. The main 
observables taken into account for composition studies at 
KASCADE-Grande are the shower size (N e ) and muon 
size (N^). For all the methods it is crucial to verify 
the sensitivity of the observables to different primary 
particles and the reproducibility of the measurements 
with the hadronic interaction model in use as a function 
of sizes and the atmospheric depth. Two examples how 
we approach the composition studies are discussed in 
the following: 

A. Muon densities 

In addition to the validity checks of the hadronic 
interaction models the reconstruction of the local muon 
density at a certain distance to the shower core gives a 
sensitivity to changes in the elemental composition 031 
(see figs. [7] and 0). The data as shown in figure [8] can 
be used to study a possible transition of the primary 
composition with increasing total number of electrons, 
but conclusive results are presently not possible due 
to the strong dependence on the hadronic interaction 
model. However, the analysis done with both models 
show no abrupt change in the composition in the entire 
energy range. EPOS 1.61 would require a very light 
abundance of primary particles in order to fit the data. 
QGS Jet-II could fit the data with an intermediate primary 
abundance. The change in slope seen in figure [8] for 
logio(N e ) < 6.0 corresponds to the threshold of the 
experiment and the fact that both data and simulation 
show the same behavior illustrates the good level of 
understanding of our detector. 

B. Electron-muon number ratio 

In the second example, the total number of electrons 
N e and the total number of muons of each recorded 
event are considered and the distribution of N^/N e is 
studied in different intervals of N e (corresponding to 
different energy intervals) and zenith angle l20l . 

KASCADE-Grande data are chosen, at first, in an 
electron size range providing full reconstruction effi- 
ciency and high statistics: 6.5 < logio(N e ) < 6.75 
in 0° < 6 < 24° (equivalent to a primary energy 
just below 10 17 eV). The same event selection is made 
on the simulated QGSJet-II data sets at disposal for 
each cosmic ray primary (p, He, C, Si, Fe) simulated 
with an energy spectrum with slope 7 = —3. The 
experimental distribution of the observable N^/N e is 
taken into account and fitted with a linear combination 
of elemental contributions from simulations. By this way 
the mean and the width of the distributions are taken into 
account. 

It was found that by using only two primary masses 
for the fit, the measured distributions can not be de- 
scribed. Fig. [321 shows the fit performed with a combi- 
nation of three elements: Protons, carbon and iron which 
is matching the light elements with the heaviest element. 
It can be seen that a three elements combination is well 
fitting the data, but it is also possible to use another 




Fig. 12: The KASCADE-Grande data in 6.5 < lg(N e ) < 
6.75 (0° < < 24°) described by Protons, Carbon and 
Iron primaries. 



— data 




Fig. 13: The KASCADE-Grande data in 7.25 < 
lg(N e ) < 7.5 (0° < < 24°) described by Protons, 
Helium and Iron primaries. 



combination of primaries, e.g. protons, helium and iron. 
However, protons and iron primaries are needed to 
describe the edges and the width of the data distribution. 
These results show that the QGS Jet model describes well 
the shape and the tails of the experimental distribution. 
Selecting the KASCADE-Grande experimental data for 
higher values of the electron size N e means to choose 
showers that were generated by higher energy events 
(fig. [13] is for energies slightly above 10 17 eV). Also in 
this case (as well as for other ranges in zenith angle), 
it is found that the model reproduces the data and 
the minimization of the N^/Nq distribution with three 
chemical components still provides a good result. 

Summarizing, by such first composition studies it is 
seen that, by exploiting a chi- square minimization of 
a linear combination of different simulated primaries, 
the KASCADE-Grande N^/Ne distributions can only be 
fitted using at least three primary elements. QGSJet- 
II, the hadronic interaction model in use, can fairly 
well reproduce the data and, in particular, the tails of 
the distributions, that represent a main constraint being 
related to the lightest and heaviest cosmic ray primaries. 
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Fig. 14: Comparison of the all-particle energy spectrum obtained with KASCADE- Grande data to results of other 
experiments. 



One has to remark, that using another hadronic inter- 
action model would of course lead to significant changes 
in the relative abundances of the elemental groups as 
different models predict different means. But, as the 
width of the distributions are more or less independent of 
the interaction model, the fact that at least three primary 
mass groups are needed to describe the data is valid. 

VII. Conclusions 

Applying various different reconstruction methods to 
the KASCADE-Grande data the obtained all-particle 
energy spectra are compared for cross-checks of recon- 
struction, for studies of systematic uncertainties and for 
testing the validity of the underlying hadronic interaction 
models. The resulting energy spectra are consistent 
and in the overlapping energy range in a very good 
agreement to the spectrum obtained by the KASCADE, 
EAS-TOP, and other experiments (fig. [14]). A wealth 
of information on individual showers is available with 
KASCADE-Grande. This makes it possible to recon- 
struct the all-particle energy spectrum with high preci- 
sion, as well as to investigate the elemental composition, 
to test the hadronic interaction models, and to study 
cosmic ray anisotropies. All these studies are under way 
and results are expected in near future. 
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